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Abstract:  A  wavelet-transform  based  strategy  is  proposed  for  the  power  management  of  hybrid  electric  vehicles  (HEV)  with  multiple 
on-board  energy  sources  and  energy  storage  systems  including  a  battery,  a  fuel  cell,  and  an  ultra-capacitor.  The  proposed  wavelet 
transform  algorithm  is  capable  of  identifying  the  high  frequency  transient  and  real  time  power  demand  of  the  HEV,  and  allocating 
power  components  with  different  frequency  contents  to  corresponding  sources  to  achieve  an  optimal  power  management  control 
algorithm.  By  using  the  wavelet  decomposition  algorithm,  a  proper  combination  can  be  achieved  with  a  properly  sized  ultra-capacitor 
dealing  with  the  chaotic  high-frequency  components  of  the  total  power  demand,  while  the  fuel  cell  and  battery  deal  with  the  low  and 
medium  frequency  power  demand.  Thus  the  system  efficiency  and  life  expectancy  can  be  greatly  extended.  Simulation  and 
experimental  results  validated  the  effectiveness  of  wavelet-transform  based  power  management  algorithm. 

Keywords:  Hybrid  electric  vehicle,  wavelet  transform,  wavelet  decomposition,  power  management,  power  demand,  power  flow. 

I.  Introduction 

Fuel  cells  (FC’s)  have  been  considered  as  an  alternative  source  of  power  for  vehicle  propulsion.  As  an  alternative  power 
source,  FC’s  are  environment  friendly  because  they  do  not  generate  any  harmful  pollution.  As  energy  converters,  FC’s  are  also 
more  efficient  than  internal  combustion  engines  (ICE)  [1-2].  However,  high  cost  and  poor  transient  performance  have  restrained 
FC’s  from  being  used  in  large-scale  vehicular  applications.  Due  to  the  unidirectional  power  flow  characteristics  of  the  FC’s,  the 
energy  from  regenerative  braking  of  a  vehicle  can  not  be  handled  by  FC’s,  hence  vehicle  efficiency  enhancement  is  not  possible 
through  these  devices  alone.  In  addition,  the  size  and  lifetime  of  an  FC  may  be  dramatically  impacted  by  the  chaotic  and  rapid 
power  demand  variations  of  an  HEV.  Therefore,  hybridization  of  fuel  cells  with  other  energy  storage  devices,  such  as 
ultra-capacitors  and  batteries,  are  often  necessary  [3-5]. 

For  systems  containing  multiple  power  sources  or  multiple  energy  storage  systems,  an  appropriate  power  management 
algorithm  is  indispensible.  Ref.  [6]  introduced  a  fuzzy-logic-supervisory-based  power  management  strategy  for  a  fuel 
cell/ultra-capacitor/battery  combined  electric  vehicle.  An  adaptive  control  method  for  power  sharing  in  hybrid  FC/battery  power 
sources  was  presented  in  [7].  These  power  management  strategies  have  been  proven  effective  in  dealing  with  system  efficiency, 
but  they  fail  to  address  the  lifetime  of  the  fuel  cell  stack  due  to  negative  impact  of  frequent  and  rapid  load  demand  variations  in 
real  world  drive  cycles. 

Ultra-capacitors  (UC’s)  offer  significantly  better  energy  densities  than  do  conventional  capacitors.  They  are  also  highly 
suitable  for  the  bulk  of  the  transient  power  demands  during  acceleration  and  braking  processes  since  the  charge/discharge  time  of 
an  UC  can  vary  in  a  wide  range  [8-10],  which  ensures  that  the  high  efficiency  and  high  durability  of  the  UC's  are  fully  employed, 
and  the  lifetime  of  the  FC  is  extended  because  the  UC’s  can  handle  lots  of  high  and  fast  power  outputs.  Also,  if  there  exists  a 
battery  pack  in  the  HEV,  the  lifetime  of  the  battery  is  extended  as  a  result  of  not  being  subjected  to  many  high-current  charges  or 
discharges.  Delay  of  battery  degradation  is  a  major  advantage  for  the  combined  battery/ultra-capacitor  energy  storage  system 
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(ESS).  The  wavelet  transform  algorithm  shows  satisfactory  performance  in  analysis  and  acquisition  of  high-frequency  transients 
[11-14]  (e.g.  transient  power  demand  variations  during  the  drive  cycles).  So  a  wavelet-decomposition  algorithm  is  introduced  in 
this  paper  for  the  power  management  of  an  HEV,  where  the  UC  deals  with  rapid-variation  components  of  the  total  power  demand 
of  a  HEV. 

The  incorporation  of  a  battery  with  high  energy  density  to  the  system  not  only  solves  the  problem  of  unidirectional 
low-frequency-component  power  flow  with  the  FC/UC  configuration,  but  also  reduces  the  size  of  the  FC  module  since  the  output 
power  is  shared  between  the  FC  and  battery.  Moreover,  the  output  voltage  of  the  battery  is  less  influenced  by  the  load  compared 
to  a  FC  [15].  Therefore  the  battery  can  help  stabilize  the  DC  bus  voltage  and  also  make  the  power  system  more  controllable. 

In  this  paper,  an  advanced  vehicle  simulator  [27]  tool  is  employed  for  modeling  the  power  demand  of  the  HEV.  The  Orange 
County  Cycle  (OCC)  load  power  profile  is  chosen  as  an  example  to  demonstrate  the  power  demand  during  the  drive  cycle.  First, 
the  OCC  power  demand  is  decomposed  into  high  and  low  frequency  components  by  utilizing  the  wavelet  transform.  Then,  the 
high-frequency  components  are  assigned  to  the  UC  bank.  The  FC  stack  and  battery  pack  share  the  low-frequency  components. 
The  battery  absorbs  the  negative  (regenerative)  portion  of  power  demand  and  the  FC  module  carries  the  major  portion  in  the 
positive  part  (including  load  power  and  battery/UC  charging).  Finally,  simulation  [28]  and  real  time  environment  [29]  based 
experimental  results  verify  that  the  power  source  modeling  is  correct,  and  the  proposed  wavelet-transform  power  management 
algorithm  is  valid. 


II.  System  Configuration 

Fig.  1  shows  a  typical  hybrid  energy  system  consisting  of  an  FC,  a  battery  and  a  UC.  A  unidirectional  DC/DC  converter  and 
two  bi-directional  DC/DC  converters  are  directly  connected  to  the  FC,  battery  and  UC,  respectively.  The  power  demand  of  the 
HEV  can  be  satisfied  by  controlling  the  direction  and  amount  of  the  output  power  (power-oriented-control)  of  each  energy 
component.  The  optional  bi-directional  DC/DC  converter  placed  between  the  DC  bus  and  the  load  is  for  the  purpose  of  voltage 
regulation  to  keep  the  DC  load  voltage  stable.  The  required  power  of  the  FC,  battery  and  UC  is  calculated  by  the  proposed 
wavelet-transform  based  algorithm. 


PWM  switch  signal  for  FC  power  tracking 


Switch  signal  (on-status  means  charging  for  Battery) 


Fig.  1.  System  configuration  of  the  multi-source  hybrid  vehicle  power  system. 
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The  terminal  voltage  and  output  current  of  each  power  source,  which  help  to  obtain  the  instantaneous  power,  are  sampled  and 
provided  as  inputs  to  the  data  acquisition  (DAQ)  and  control  system  for  comparison  with  the  desired  power.  The  DC  voltage  on 
the  load  side  is  sampled  and  fed  back  to  the  DC  load  voltage  regulator  located  in  the  DAQ  and  control  system.  The  PWM  signals 
generated  from  the  PID  controllers  are  used  to  control  the  four  DC/DC  converters  respectively  for  power  tracking  and  voltage 
regulation.  If  there  are  only  three  DC/DC  converters  without  the  optional  DC/DC  converter  connected  between  the  load  and  DC 
bus,  the  DC  load  voltage  can  be  stabilized  by  choosing  the  voltage-oriented-control  for  the  converter  connected  to  the  battery 
considering  its  voltage  robustness,  and  power-oriented-control  for  the  other  two  converters  connected  to  the  FC  and  UC.  If  the 
state  of  charge  (SOC)  of  the  battery  or  UC  falls  below  a  predetermined  threshold,  the  switch  between  the  FC  and  the  battery/UC 
is  closed  until  the  voltage  returns  within  the  desired  limits. 

The  output  current  from  the  FC  module  is  limited  by  its  maximum  value.  The  battery  and  UC  are  initially  charged  to  sufficient 
voltage  to  ensure  that  the  voltage  difference  between  the  FC  and  the  battery  or  the  UC  is  within  an  acceptable  range  when  the 
switch  turns  on.  This  helps  ensure  that  the  surge  current  is  within  allowable  limits  and  the  membrane  of  the  FC  module  is  not 
damaged.  The  specific  control  process  will  be  discussed  in  detail  in  Section  IV-C. 

III.  Modeling  of  Power  Sources 

In  order  to  develop  and  validate  the  proposed  power  management  algorithm,  the  models  of  each  component  need  to  be 
developed. 

A.  Modeling  of  Lithium-Ion  Battery 

Lithium-ion  batteries  are  idea  candidate  for  use  in  hybrid  electric  vehicles  due  to  their  high  energy  density,  high  power  density, 
and  long  lifetime  [16-18].  This  paper  uses  a  lithium-ion  battery  as  a  power  source  capable  of  absorbing  the  slow-variation 
stationary  regenerative  energy  and  assisting  the  FC  to  provide  adequate  power  to  the  load. 


V 

Fig.  2.  Equivalent  circuit  model  of  the  lithium-ion  battery. 

The  equivalent  electrical  schematic  of  a  lithium-ion  battery  is  shown  in  Fig.  2.  It  consists  of  an  equilibrium  potential  E  ,  an 
internal  resistance  having  two  components  R]  and  R2  and  an  effective  capacitance  C  characterizing  the  transient  response 
of  charge  double  layers. 

The  equilibrium  potential  of  the  battery  depends  on  the  temperature  and  the  state  of  discharge  (SOD).  The  reference  discharge 
curve,  reference  rate  and  reference  temperature  are  chosen  based  on  experimental  data.  The  equilibrium  potential  as  a  function  of 
SOD  is  expressed  as  an  nth-order  polynomial.  A  potential  correction  term  is  used  to  compensate  for  the  variation  resulting  from 
the  ambient  temperature  change.  Thus,  the  following  equations  are  obtained  to  express  the  relationship  among  the  potential, 
terminal  voltage  and  SOD  [19]. 

E[i(t),  T  (0,  t ]  =  v[i(t),  T  (0,  t  ]  -  RJ(t)  (1) 

V[i(t),  T(t),  fJ  =  X  q  SODk  [i(t),  T(t),t]  +  A  E(T) 

k=  0 


(2) 
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SOD[i(t),  T{t),  t]  =  —  r  oc[i(t)\  ■  j8[T(t)]  ■  i(t)dt 

C  Jo 


(3) 


where  Ck  is  the  coefficient  of  the  kth  order  term  in  the  polynomial  expression  for  the  reference  curve  and  Cr  is  the  battery 

capacity;  A E  expresses  the  potential  correction  term. 

The  battery  temperature  change  is  determined  by  the  thermal  energy  balance  expressed  as  [20] 

m.c  ■<i^  =  i(t)2  ■Rl+^[v(t)-E[i(t),T(t),t]-i(t)Rlf  -hcA[T(t)-Ta]  (4) 

p  dt  R2 

where  m  is  the  battery  weight;  A  is  the  surface  area;  Cp  is  heat  capacity;  It  is  cooling  coefficient;  T  represents  the 

ambient  temperature. 

B.  Modeling  of  Fuel  Cell 

Proton  exchange  membrane  fuel  cells  (PEMFC’s)  are  considered  to  be  promising  for  application  as  distributed  generation 
sources  due  to  high  efficiency  and  compactness  [1-2],  [21-24].  The  PEMFC  is  also  a  primary  candidate  for  environment-friendly 
vehicles. 

The  nomenclatures  of  the  PEMFC  are  as  follows: 

B  ,  C  :  Constants  to  calculate  the  activation  over  voltage  [A'1]  and  [V];  E  ^  :  Nernst  cell  voltage  [V];  -Edcell  Overall  effect 
of  the  fuel  and  oxidant  delay  [V];  EQce]1  Open  cell  voltage  [V];  ^  Standard  reference  potential  per  cell  [V];  F 

Faraday's  constant  [C/kmol];  I  FC  current  [A];  kE  Empirical  constant  in  calculating  open  cell  voltage  [V/K];  kH 
Hydrogen  valve  molar  constant  [kmol/(atm-s)];  kQ  Oxygen  valve  molar  constant  [kmol/(atm-s)];  kr  Modeling  constant 


[kmoF(s-A)];  N0  Number  of  FCs  in  series;  pH  Hydrogen  partial  pressure  [atm] ;  p()  Oxygen  partial  pressure  [atm] ; 
Hydrogen  input  flow  [kmol/s];  <7^q  Required  hydrogen  flow  for  load  change  [kmol/s] ;  Oxygen  input  flow  [kmol/s]; 
tp)1'  Required  oxygen  flow  for  load  change  [kmol/s];  R  Gas  constant  [J/(mol  K)];  rH  0  Hydrogen-oxygen  flow  ratio;  Rohm 
Internal  resistance  [Q];  T  FC  stack  temperature  [K];  U  Utilization  rate;  Vact  Activation  over  voltage  [V];  V|)hrn  Ohmic 


over  voltage  [V];  Vout  FC  output  voltage  [V];  A  Constant  factor  in  calculating  Edcell  [Q];  Te  Overall  flow  delay  [s]; 


Th  Hydrogen  time  constant  [s];  T()  Oxygen  time  constant  [s]. 


At  the  anode,  three  factors  (hydrogen  input  flow,  hydrogen  output  flow  and  hydrogen  flow  that  reacts)  determine  the  effective 
hydrogen  partial  pressure.  The  final  expression  for  the  hydrogen  partial  pressure  is  described  as  [22] 

1  IK 


Ph, 


1  +  T„  5 


«  -2 krI) 


(5) 
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Similarly,  the  expression  for  the  oxygen  partial  pressure  is  obtained  as 

l/*o, 

Po,  = 


1 +To,S  2 


«  -KD 


(6) 


The  relationship  among  the  FC  output  voltage  V.  ,  Nernst  cell  voltage  EccH ,  activation  voltage  V!  =51n(C7)  and 
overall  ohmic  voltage  V()hm  =  IRohm  is  delineated  as  [23]: 


where 


V^N^-BMCD-IR 


ohm 


^cell  ^O.cell  ' 


RT 


kE(T-29S)  +  —  \n[pH2-(p02r]-Ei 


d.cell 


(7) 

(8) 


^d.cell  (s)  =  AeI(s)^~  (9) 

Tes+l 

It  should  be  noticed  that  the  fuel  cell  voltage  has  to  be  kept  at  or  below  0.9V  otherwise  the  Pt  catalyst  oxidizes  and  the 
platinum  migrates.  In  this  paper,  the  no-load  voltage  of  the  studied  fuel  cell  is  0.6V.  That  is  to  say,  at  any  load,  the  fuel  cell 
voltage  will  be  at  or  below  0.6V,  so  the  aforementioned  requirement  is  satisfied.  Specific  parameters  and  numerical  analysis  will 
be  shown  in  Part  A  of  Sec.  V. 

A  feedback  from  the  FC  module  current  is  applied  to  control  the  hydrogen  and  oxygen  flow  rates  according  to  the  power 
demand 


req 

^  2  FU 


rH-0 


N0I 

Wo  FU 


The  block  diagram  of  the  PEMFC  model  based  on  the  above  equations  is  shown  in  Fig.  3. 


i 


(10) 


Fig.  3.  PEMFC  model. 
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Fig.  4.  Equivalent  circuit  model  of  the  ultra-capacitor. 

The  classical  equivalent  electric  model  [8,9]  of  a  UC  cell  is  introduced  as  shown  in  Fig.  4.  The  circuit  parameters  are  shown  as 
follows: 

C  :  Capacitance  [F];  Cucbank  :  Total  capacitance  of  the  UC  bank  [F];  E  :  Amount  of  energy  flowed  from  the  UC  [W-s];  Ylx : 
Number  of  capacitors  in  series;  n2 :  Number  of  UC  branches  in  parallel;  Repr-  Equivalent  parallel  resistance  | Q  | ;  Resr  : 
Equivalent  series  resistance  [Q];  Vt) :  Initial  voltage  [V];  Vf  :  Final  voltage  [V]. 

As  depicted  in  Fig.  3,  the  series  resistance  Resr  is  the  charging/discharging  resistance,  and  the  parallel  resistance  Repr 


represents  the  self-discharging  losses.  However,  Repr  only  has  impact  on  the  long-term  energy  storage  performance. 


The  state  of  charge  (SOC)  of  the  UC  can  be  described  as  a  percentage  of  the  rated  energy  capacity,  which  depends  on  the 
terminal  output  voltage.  The  energy  flowing  out  from  the  UC  is  directly  determined  by  the  capacitance  and  the  voltage  change  as 
expressed  in  ( 1 1 )  [9,10], 

E  =  ~C(Vq  -Vf)  (11) 

2 

The  real  UC  bank  can  be  modeled  by  using  multiple  UC  cells  in  parallel  and  series.  The  total  resistance  and  capacitance  of  the 
UC  bank  are  given  by  [9]: 


R  —  // 

ivUCbank  "l 


IV.  Wavelet-Transform  Based  Power  Management  Strategy 
There  have  been  a  variety  of  strategies  capable  of  managing  the  power  among  different  power  sources.  Fuzzy  logic  is  well 
known  to  enhance  conventional  system  design  with  engineering  expertise,  and  the  use  of  fuzzy  logic  can  help  to  circumvent  the 
requirement  for  rigorous  mathematical  modeling.  Fuzzy  logic  supervisory  control  applied  in  multiple  power  sources  of  a  HEV 
can  avoid  the  low  efficiency  operation  region  of  every  power  source  and  achieve  high  efficiency  without  degrading  the 
mechanism  performance.  Ref.  [6]  proposed  a  power  management  strategy  based  on  the  fuzzy  logic  supervisory  system  for  a  fuel 
cell/ultra-capacitor/battery  combined  electric  vehicle,  whose  objectives  are  to  improve  the  performance  of  individual  power 
sources  and  to  ensure  high  power  quality  desirable  for  a  good  performance  for  electrical  loads.  Another  power  management 
strategy,  adaptive  control  method,  was  presented  in  [7]  for  hybrid  FC/battery  power  sources.  In  this  adaptive  control  system,  the 
hybrid  power  sources  are  optimized  so  that  the  fuel  cell  provides  the  energy  required  by  the  load,  while  the  battery  satisfies  the 
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peak  power  demands  of  the  load.  The  output  current  of  the  fuel  cell  is  adjusted  according  to  the  charge  level  of  the  battery.  By 
this  means,  the  fuel  cell  can  achieve  high  efficiency,  and  the  battery  will  not  run  down  to  depletion. 

The  power  management  strategies  mentioned  above  have  been  proven  effective  in  dealing  with  system  efficiency.  However, 
they  fail  to  address  the  lifetime  of  the  fuel  cell  stack  and  battery  due  to  the  dramatically  negative  impact  of  the  frequent  and  rapid 
load  demand  variations  in  real  world  drive  cycles.  The  lifetime  of  the  fuel  cell  is  a  critical  factor  for  HEV  economy.  So  the 
questions  of  how  to  avoid  the  rapid-variation  and  large-differential  power  demand  for  the  fuel  cell  and  how  to  appropriately 
distribute  the  power  demands  to  power  sources  according  to  their  natural  characteristics  arise.  The  wavelet  transform  (WT)  can 
provide  a  potential  solution.  The  WT  decomposes  an  original  signal  into  components  at  different  positions  and  scales.  It  can  be 
applied  to  extract  signal  information  in  both  time  and  frequency  domains.  Wavelets  with  the  property  of  localization  are  suitable 
for  the  analysis  of  non-stationary  and  transient  signals  [11-14].  Therefore,  WT  can  be  employed  as  a  very  effective  tool  for  the 
decomposition  and  detection  of  the  transient  features  in  1-D  seismic  data,  and  similar  data  such  as  the  power  demand  of  the 
vehicle  in  this  study.  We  extract  the  power  transients  from  the  total  power  demand  and  deliver  them  to  the  UC  since  the 
charging/discharging  time  of  the  UC  can  vary  from  several  milliseconds  to  minutes.  This  measure  protects  the  FC  from  the 
membrane  damage  resulting  from  the  quick  power  demand  variations.  On  the  other  hand,  the  battery  absorbs  the  negative  part  of 
the  rest  power  demand  and  shares  the  positive  part  with  the  FC.  The  FC  does  not  have  to  provide  the  high  peak  power 
occasionally  demanded  by  the  load.  As  a  result  of  using  this  WT-based  power  management  strategy,  the  high  efficiency  of  the  FC 
can  be  achieved  and  the  lifetime  of  the  FC  can  be  extended. 

A.  Algorithm  Analysis 

The  parameters  used  for  WT-based  algorithm  analysis  in  this  section  are  as  follows: 

G  Decomposition  transfer  function  matrix;  H  Reconstruction  transfer  function  matrix;  K  Constant  determined  by  Iff  \  t 


Time;  U  Position  factor;  W  Wavelet  coefficients;  X  Original  signal;  X  Scale  factor;  If/  Mother  function 
The  Morlet-Grossmann  definition  of  continuous  wavelet  transform  (CWT)  for  a  1-D  signal  x(t)  is  described  as 

W(X,u)=\  — — )d t,  X>0,  ueR 

Jr  sJX  X 


(13) 


where  X  is  the  scale  parameter  and  U  is  the  position  factor.  \j/(l )  is  a  mother  function  in  both  time  domain  and  frequency 


domain  called  the  mother  wavelet.  This  mother  wavelet  is  for  the  purpose  of  providing  a  source  function  to  generate  the  daughter 
wavelets  which  are  simply  the  translated  and  scaled  versions  of  the  mother  wavelet.  The  inverse  continuous  wavelet  transform 
can  be  exploited  to  recover  the  original  signal  x(t), 


x(t)=^-LLwa-u)^( 


u  dXd  u 


4xYK  x  '  x 


(14) 


As  the  most  popular  mother  wavelet,  the  Haar  wavelet  has  the  shortest  filter  length  in  the  time  domain  compared  to  other 
wavelets.  Also,  the  Haar  wavelet  is  the  simplest  possible  wavelet  with  a  nice  feature  that  the  WT  is  equal  to  its  inverse  [13].  With 
these  features,  the  Haar  wavelet  makes  the  decomposition  calculation  much  simpler  than  other  wavelets  do  for  realization  of  the 
wavelet  transform  strategy  in  a  realistic  system.  This  can  simplify  the  program  composition  and  further  improve  the  code 
execution  efficiency.  Meanwhile,  the  function  of  extracting  transients  can  still  be  implemented  well  without  degradation. 
Consequently  the  Haar  WT  is  used  in  this  paper  for  load  profile  decomposition  to  obtain  the  transients  that  are  assigned  for  the 
ultra-capacitor.  The  Haar  wavelet  is  expressed  as 
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Unlike  CWT,  the  discrete  wavelet  transform  (DWT)  is  used  to  decompose  a  discretized  signal  into  different  resolution  levels. 
The  DWT  and  its  inverse  are  given  by  the  following  two  equations  respectively. 


W(A,u)={  x(t) 

JR 


t-u 

~Y~ 


)d  t 


A  =  2j ,  u  =k 2j ,  j,ke  Z 


x(t)  =  ZYjW^k)VLk(t) 

jeZ  keZ 


(16) 

(17) 


As  shown  in  Fig.  5,  by  using  the  Haar-wavelet  based  A'-channel  filter  bank  [  14],  the  original  signal  x(t)  can  be  decomposed  by 
a  low-pass  filter  H0(z )  and  high-pass  filters  H k  (z) ,  k  =  1, 2, •  •  •  N  —  1 ,  into  a  reference  signal  and  other  detailed  signals 

respectively.  The  reconstruction  filters  Gk(z)  ,k  =  0, 1, •  •  •  iV  —  1  are  placed  for  signal  reconstruction. 


Fig.  5.  Structure  of  Haar- wavelet-based  iV-channel  filter  bank. 
The  decomposition  transfer  function  matrix  is  expressed  as 
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where  z  is  the  operator  of  the  Z-transform. 
The  reconstruction  filters  are  given  by  [14] 


(18) 
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The  number  of  wavelet  decomposition  levels  has  to  be  considered  properly  complying  with  calculation  simplicity  and 
frequency  restrictions  as  well  in  a  realistic  system.  This  number  is  analyzed  and  eventually  determined  to  be  three.  Detailed 
analysis  can  be  found  in  Part  C  of  this  section.  Consider  N=2,  then  the  decomposition  and  reconstruction  filters  become 

[//j(z),  H0(z)]T  and  [Gj(z),  G0(z) \'  ,  respectively.  The  three-level  Haar  wavelet  decomposition  and  reconstruction  are 

used  for  input  signal  x(t).  As  shown  in  Fig.  6,  the  down-sampling  and  up-sampling  methods  are  employed  in  the  decomposition 
and  reconstruction  processes  respectively.  The  data  size  reduces  by  half  in  down-sampling  operations  while  it  doubles  in 
up-sampling  operations. 


Fig.  6.  Three-level  Haar  decomposition  and  reconstruction  diagram. 

B.  Power  Demand  Sharing  Structure 

Advanced  simulation  tool  is  employed  to  help  set  up  a  power  demand  model  of  the  proposed  power  management  system  in  an 
HEV.  The  Orange  County  Cycle  (OCC)  is  chosen  as  a  desired  cycle  for  modeling.  The  speed  curve  during  the  OCC  is  drawn  in 
Fig.  7.  The  power  demand  of  the  OCC  driving  cycle  is  shown  in  Fig.  8.  For  the  convenience  of  design,  the  load  profile  is 
considered  as  the  power  demand.  It  can  be  seen  from  Fig.  8  that  the  load  profile  consists  of  many  transients  derived  from  the 
speed  variations.  By  using  the  Haar-wavelet-transform  based  scheme,  this  power  demand  is  divided  into  three  parts  for  power 
distribution,  which  can  take  advantage  of  the  three  sources  (i.e.  the  FC,  UC  and  battery).  The  WT  can  properly  deal  with  the 
nonlinear  and  uncertain  transients. 
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Fig.  7.  Speed  curve  of  OCC. 


Fig. 8.  Load  profile  of  OCC. 


The  principle  depicted  in  Fig.  6  is  considered  as  a  basis  for  the  power  sharing  strategy.  H0(z)  is  a  low-pass  filter  and 


tfj(z)  is  a  high-pass  filter.  The  original  signal  X (n)  equals  the  predetermined  power  demand  of  the  OCC.  As  shown  in  Fig.  6, 
after  three-level  low-pass  processing,  the  smooth  slow-variation  signal  X0(n)  is  obtained  as  the  power  demand  for  the  FC  and 


battery.  Other  transients  ( Xl (n)  +  x2 (n )  +  A',  (n) )  are  distributed  to  the  UC  as  mentioned  above.  However,  X0 (n)  contains  a 


negative  part,  which  the  FC  cannot  accept.  Thus,  the  battery  absorbs  the  slow-variation  power  with  the  direction  from  the  load  to 
power  sources.  The  battery  also  helps  absorb  a  portion  of  the  positive  low-frequency  components,  while  the  FC  deals  with  the 
significant  portion.  So  the  power  sharing  method  is  expressed  as  below: 

["0.6  x0(n) 

0 


p  =  . 

'  fc 


if  x0  (n)  >  0 
otherwise 


(20) 


[  Battery 


\0AxQ{n) 
I  x0(n) 


if  x0(n)>  0 
otherwise 


(21) 


1  uc 


:  xx(n)  +  x2(n)  +  x3(n) 


(22) 


Based  on  this  sharing  strategy,  the  power  demands  for  the  FC,  battery  and  UC  are  obtained  and  shown  in  Fig.  9  (a),  (b)  and  (c), 
respectively.  It  can  be  observed  that  the  power  demand  curve  for  the  FC  has  global  low  frequency  using  the  Haar  wavelet 
decomposition,  while  at  some  time  instants,  fast  power  demand  for  the  FC  still  exists.  However,  the  number  of  large  and  fast 
power  variations  reduces  dramatically,  and  the  required  surge  power  for  the  FC  also  decreases  largely  because  of  the  assistance 
of  the  battery  (compare  Fig.  9(a)  with  Fig.  8).  Consequently,  the  FC  connected  to  a  power  converter  with  appropriate  capacity 
and  fast  response  can  handle  this  power  demand.  The  lifetime  of  the  FC  is  extended  consequently,  which  is  a  major  advantage  for 
utilization  of  the  Haar  wavelet  transform  in  HEV  power  management. 
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(a)  Power  demand  for  FC. 


(b)  Power  demand  for  battery. 


(c)  Power  demand  for  UC. 


Fig.  9.  Power  components  for  the  three  power  sources. 


C.  Data  Acquisition  and  Control  System 

For  the  three-level  wavelet  transform  (see  Fig.  6),  the  power  variation  frequency  of  the  FC  and  battery  is  only  1/8  of  that  of  the 
UC,  so  only  eight  sampling  power  points  need  to  be  foreknown  for  every  calculation  in  power  demand  decomposition.  The 
power  demand  studied  in  this  paper  is  already  known  due  to  the  predefined  drive  cycle.  By  using  the  power  prediction  algorithm 
[25]  considering  the  vehicle  speed  and  motor  operation  period,  we  can  easily  implement  the  real-time  power  management  for  a 
running  HEV.  The  calculation  size  is  not  large  enough  so  that  even  a  single-chip-machine  can  handle  it. 

As  shown  in  Fig.  10,  an  isolated  gate  bipolar  transistor  (IGBT)  SI  is  employed  in  the  topology  of  the  unidirectional  boost 
DC/DC  converter  for  power  tracking  of  the  FC  as  mentioned  in  Sec.  II.  The  power  from  the  fuel  cell  to  the  DC  bus  can  be 
determined  by  inputting  proper  PWM  drive  signals  to  S 1 .  By  adding  IGBT  S2  and  diode  D 1 ,  the  unidirectional  DC/DC  converter 
changes  to  a  bi-directional  converter,  as  shown  in  Fig.  1 1 .  For  power  tracking  applications  of  the  battery  and  UC  as  described  in 
Sec.  II,  the  battery  and  UC  are  located  on  the  Viow  side,  and  the  other  side  Vhigh  is  connected  to  the  DC  bus.  So  the  voltage  from 
the  power  sources  is  boosted  on  the  DC  bus  with  the  aim  of  voltage  regulation,  the  DC  bus  is  placed  on  the  Vhigh  side  and  the 
load  is  linked  to  the  Viow  side.  The  DC  load  voltage  is  regulated  to  a  constant  approaching  the  power  sources'  voltage  values. 


D2 


L 

- N - 
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Fig.  10.  Unidirectional  DC/DC  converter  connected  between  FC  and  DC  bus. 
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Fig.  11.  Bi-directional  DC/DC  converter  topology. 
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The  presence  of  DC/DC  converters  may  increase  the  system  cost  and  volume  compared  to  the  direct  connected  hybrid  systems 
[26].  But  by  using  DC/DC  converters,  the  FC  module  size  can  be  optimized.  The  FC  is  also  protected  from  high-frequency  and 
surge  load  power  demand.  Therefore  the  life  expectancy  can  be  extended.  Meanwhile,  the  regenerated  energy  is  properly  restored 
due  to  the  existence  of  UC  and  battery. 

In  this  proposed  power  management  system  (see  Fig.  1),  the  three  power  sources  (i.e.  FC,  battery  and  UC)  and  corresponding 
power  flow  control  are  independent.  The  required  power  demand  suitable  for  the  characteristics  of  an  individual  power  source 
can  be  met  precisely.  The  specific  control  concepts  are  described  below: 

•  The  three-level  Haar  wavelet  transform  is  utilized  for  decomposition  of  the  HEV  power  demand. 

•  The  PEMFC  module  deals  with  the  majority  of  the  positive  portion  of  the  low-frequency  components  of  power 
demand  derived  from  the  wavelet  decomposition.  The  power  flow  from  the  FC  is  determined  by  switch  duty  cycle  of 
the  unidirectional  DC/DC  converter.  The  duty  cycle  is  obtained  by  the  DAQ  and  control  system. 

•  The  battery  satisfies  the  negative  portion  (during  regeneration)  and  some  of  the  positive  portion  of  the  low-frequency 
components  of  power  demand  derived  from  the  wavelet  decomposition.  The  battery  not  only  absorbs  the 
slow-variation  power  from  the  load  but  also  helps  the  FC  provide  the  required  power  to  reduce  the  burden  of  the  FC. 

•  All  the  high-frequency  components  of  the  power  demand  are  met  by  the  ultra-capacitor  (UC).  The  amount  and 
direction  of  required  fast-variation  power  is  determined  by  properly  controlling  an  individual  switch  existing  in  the 
bi-directional  DC/DC  converter.  The  switch  status  and  switch  duty  cycle  are  outputs  of  the  DAQ  and  control  system. 

•  When  the  SOC  of  the  battery  drops  to  0.5  because  of  over-discharging,  switch  Sbat  in  Fig.  1  turns  on  (i.e.  connected  to 
the  upper  diode  attached  to  FC)  until  the  battery  voltage  reaches  the  FC  module  terminal  voltage.  During  this  process, 
the  battery  is  charged  through  the  diode  Dbat  and  switch  Sbat.  At  this  time  switch  Sbat  is  turned  on  and  the  battery  is 
then  charged  from  the  FC  module.  When  the  SOC  of  the  UC  decreases  to  0.25,  switch  Sue  in  Fig.  1  turns  on  (i.e. 
connects  to  Duc)  until  the  UC  voltage  reaches  the  FC  module  terminal  voltage.  The  UC  is  then  charged  from  the  FC 
module,  thus  avoiding  the  voltage  collapse. 

•  If  the  SOC  of  the  battery  or  UC  increases  to  0.95,  the  other  two  power  sources  (the  UC  &  FC,  or  the  battery  &  FC) 
are  removed  from  the  power  train,  and  the  battery  or  UC  satisfies  all  the  power  demand  as  the  only  power  source 
until  its  SOC  drops  to  0.85. 

•  The  output  current  from  the  FC  module  is  limited  to  its  maximum  value.  Meanwhile,  the  battery  and  UC  are  initially 
charged  to  sufficient  voltage,  and  the  battery  and  UC  are  also  charged  from  the  load  during  the  energy  regeneration 
process.  Thus,  the  SOC  of  the  UC/battery  will  not  decrease  widely  and  the  voltage  difference  between  the  FC  and 
battery/UC  will  be  within  an  acceptable  range  when  the  charging  switch  (Sbator  Sue)  turns  on.  These  measures  and 
conditions  guarantee  that  the  destructive  surge  current  can  be  avoided. 

V.  Simulation  Results 

A.  Verification  of  Power  Source  Models 

The  fuel  cell  parameters  are  derived  from  the  data  obtained  of  a  fuel  cell  module  [28].  According  to  the  parameters  shown  in 
Table  I,  the  open  cell  voltage  EQceU  =  u  —  kE(T  —  298)  ~  0.6V .  The  no-load  terminal  voltage  of  the  FC  module  is 

reconstructed  to  be  212V,  so  the  number  of  FC’s  in  series  is  N0  =  353  .  To  completely  satisfy  the  load  demand  shown  in  Fig. 

8(a),  two  FC  modules  are  connected  in  parallel.  The  rated  output  power  per  cell  is  27. 3W.  As  a  result,  the  rated  power  of  the  FC 
system  is  19.3kW.  Consequently,  the  required  power  in  Fig.  9(a)  can  be  completely  satisfied. 

The  lithium-ion  battery  parameters  are  shown  in  Table  II,  which  are  based  on  experimental  data  from  a  18650  lithium  battery 


Page  13  of  20 


cell.  The  nominal  voltage  of  this  battery  is  3.7V.  Accordingly  the  number  of  batteries  in  series  is  calculated  to  be  58  since  the 
terminal  voltage  of  the  battery  stack  is  212V  in  our  design.  To  obtain  the  sufficient  current  outputting  ability,  we  place  forty 
18650  lithium  battery  cells  in  parallel.  Additionally  the  nominal  capacity  of  an  18650  lithium  battery  cell  is  1.4Ah,  and  the 
reference  current  for  battery  modeling  is  set  to  1.4A.  So  the  nominal  power  of  the  battery  stack  is  12kW.  Meanwhile  the 
maximum  current  of  the  US  1 8650  battery  can  reach  a  few  tens  of  amperes.  Thus  the  maximum  instantaneous  power  of  the  stack 
exceeds  lOOkW.  The  analysis  shows  that  this  simulated  battery  stack  can  completely  meet  the  requirements  delineated  in  Fig. 
9(b). 


TABLE  I  PEMFC  MODULE  PARAMETERS 


FC  stack  temperature,  T  [K] 

338 

Number  of  FCs  in  series,  N 

353 

Faraday’s  constant,  F  [C/kmol] 

96484600 

Hydrogen  time  constant,  rH  [s] 

3.37 

Oxygen  time  constant,  th  [s] 

6.74 

Hydrogen  valve  molar  constant,  kH  fkmol/(atm-s)] 

4.22xl0"5 

Oxygen  valve  molar  constant,  kQ  [kmol/(atm  s)] 

2.11xl0'5 

Hydrogen-oxygen  flow  ratio,  rH0 

1.1679 

Standard  reference  potential  per  cell,  £0°ce|1  [V] 

0.6 

Empirical  constant  in  calculating  open  cell  voltage, 

K  rv/Ki 

0.00085 

Overall  flow  delay,  r  [s] ; 

80 

Constant  factor  in  calculating  £  ,  A  [0| 

0.00333 

Gas  constant,  R  [J/(mol  K)] 

8314.47 

Constant  to  calculate  the  activation  over  voltage,  B 

[A'1] 

0.0478 

Constant  to  calculate  the  activation  over  voltage,  C 

[V] 

0.0136 

Utilization  rate,  U 

0.85 

Internal  resistance,  R ,  [Q] 

ohm  L  J 

0.004 

The  ultra  capacitor  specifications  are  given  in  Table  III.  The  maximum  float  voltage  is  2.3V.  Therefore  105  UC  cells  are 
chosen  so  that  the  terminal  voltage  of  the  UC  bank  is  equal  to  212V.  The  total  UC  internal  resistance 

7?ucbank  =  0.00085x105  =  0.089212  and  the  total  capacitance  CUChank  =  2700/105  =  25.71F  .  The  rated  current  of  the 

UC  bank  is  400A.  Therefore  the  maximum  output  power  is  larger  than  400x212  =  84.8kW  .  This  can  adequately  satisfy  the 
power  demand  depicted  in  Fig.  9(c). 
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TABLE  II  Lithium-Ion  Battery  Parameters 


Ambient  temperature 


Parameters  referring  to  the  thermal  model 

A  =4.3E-3m2, 

m  =4 1  g. 

cp  =925J/kg/K3V/°C, 

h  t  =10W/nr/K 

Parameters  referring  to  the  equivalent  electrical  circuit 

Rt  =80mTl, 

C  =4F, 

Rn  =40mD, 

v|„  =4.2V 

Parameters  for  calculating  the  state  of  discharge 

etc  =0.93, 

charg 

a  h  =1.05, 

discharg 

0=1 

Potential  correction 

AE  =0V 

TABLE  III  Ultra-Capacitor  Parameters 

Dimensions  [mm] 

164x62x62 

Capacitance  C  [F] 

2700 

Equivalent  series  resistance  Resr  [O] 

0.00085 

Maximum  float  voltage  [V] 

2.3 

Fig.  12.  Comparison  of  the  simulated  and  measured  V-I 
characteristics  curves  of  the  FC  power  module. 


Fig.  13.  Comparison  of  simulated  and  measured 
discharging/charging  voltage  variations  of  a  lithium-ion  battery. 
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Fig.  14.  Comparison  of  simulated  and  measured  UC  voltage  during  a  period  including  charging  &  discharging. 

The  measured  and  simulated  V-I  characteristics  curves  of  the  FC  power  module,  the  charging/discharging  V-t  curves  of  the 
lithium-ion  battery  cell  and  the  UC  charge/discharge  voltage  curves  are  shown  in  Fig.  12,  Fig.  13  and  Fig.  14,  respectively.  The 
observation  that  the  simulation  results  are  very  close  to  the  measurement  data  indicates  that  the  power  source  simulation  models 
are  valid  and  accurate. 

B.  Validation  of  Power  Sharing  Algorithm 

Figs.  15-18  show  the  voltage  variations  of  the  FC  system,  battery  stack,  UC  and  DC  load  side.  Due  to  the  required  power,  the 
voltages  of  the  three  power  sources  change  accordingly.  The  battery  shows  good  voltage  stabilization  performance  in  the  whole 
drive  cycle  except  a  few  moments.  The  increase  in  assigned  power  decreases  the  FC  terminal  voltage.  The  UC  terminal  voltage 
increases  from  197V  to  215V  as  a  whole  because  of  regenerative  energy.  The  DC  voltage  on  the  load  side  remains  at  200V  with 
the  help  of  the  bi-directional  DC/DC  converter  for  voltage  regulation. 


Fig.  15.  FC  terminal  voltage. 


Fig.  16.  Battery  terminal  voltage. 
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Fig.  17.  UC  terminal  voltage. 


Fig.  18.  Regulated  DC  load  voltage. 


The  power  tracking  performance  of  the  three  power  sources  (i.e.  the  FC,  battery  and  UC)  is  shown  in  Fig.  19-21.  By 
comparing  these  power  tracking  curves  to  the  required  power  curves  shown  in  Fig.  9,  it  is  found  that  the  power  demands  are 
adequately  satisfied.  The  power  variations  are  within  the  maximum  limitations  of  the  power  sources  as  mentioned  above. 
According  to  Fig.  21,  the  UC  absorbs  the  entire  llOkJ  energy  for  a  small-scale  voltage  increase,  which  conforms  to  equation  (11). 
The  SOC  of  the  UC  during  such  a  drive  cycle  is  within  the  desired  range  (<0.95),  so  no  extra  action  is  needed.  However,  if  the 
SOC  of  the  UC  increases  to  0.95,  the  other  two  power  sources  (the  battery  &  FC)  will  be  removed  in  connection  from  the  power 
train  until  the  SOC  of  the  UC  drops  to  0.85.  The  fact  that  the  DC  output  power  trajectory  matches  the  load  profile  as  shown  in 
Fig.  8  proves  the  validity  and  correctness  of  the  proposed  power-sharing  algorithm  as  evidenced  in  Fig.  22. 
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Fig.  19.  FC  output  power. 


Fig.  20.  Charging  and  discharging  power  of  battery. 
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Fig.  21.  Charging  and  discharging  power  of  ultra-capacitor.  Fig.  22.  Total  output  power  for  load  demand. 

VI.  Experimental  Results 

A  scaled-down  experimental  setup  was  developed  in  the  laboratory  to  validate  the  proposed  algorithm.  The  power  sources 
consist  of  a  1.2kW  fuel  cell  module  (the  rated  operating  voltage  is  26V),  a  lithium-ion  battery  pack  (Vnomill;l|=25.2V,  C=32.4Ah) 
and  an  ultra-capacitor  (the  capacitance  is  95F).  For  the  sake  of  simplicity  of  experimental  verification,  a  certain  time  interval 
(97 1st  -1000th  second)  is  extracted  from  the  OCC  drive  cycle  to  be  analyzed.  Meanwhile  the  power  demand  is  scaled  down 
because  the  power  source  voltage  is  low  and  the  current  cannot  exceed  some  value  due  to  the  device  limitation.  The  power 
demand  in  the  experiment  is  1/20  of  that  in  the  real  OCC  drive  cycle,  which  is  implemented  by  an  electronic  load  (Pmax=2.6kW) 
operating  as  a  load  emulator.  Three  DC/DC  converters  controlled  by  a  real  time  controller  board  are  connected  to  three  power 
sources  respectively,  switching  at  10kHz.  The  control  program  is  compiled  and  downloaded  to  the  real  time  controller  board  for 
sending  and  receiving  signals.  Currents  and  voltages  are  monitored  and  fed  back  into  the  controller  board  where  the  power  is 
calculated  as  a  product  of  the  current  and  voltage.  The  real  instantaneous  power  (in  form  of  analog  voltage)  can  be  measured  via 
DAC  ports  on  the  real  time  controller  board. 

The  ideal  power  demand  and  the  real  power  consumed  by  the  electronic  load  are  shown  and  compared  in  Fig.  23.  The 
electronic  load  shows  good  power  tracking  precision  to  guarantee  persuasion  of  experiments.  The  designed  lkW  hardware 
prototype  is  capable  of  meeting  the  requirement  since  the  maximum  power  demand  is  600W.  The  power  demands  allocated  to 
the  fuel  cell,  battery  and  ultra-capacitor  are  calculated  using  the  proposed  wavelet  transform  and  shown  in  Fig.  24.  It  can  be  seen 
that  the  high-frequency  components  are  assigned  for  the  UC,  while  the  FC  and  battery  cooperate  with  each  other  to  share  the 
low-frequency  components  obeying  the  aforementioned  principle  (see  Eqs.  20  and  21).  Fig.  25  shows  the  measured  waveforms 
of  power  possessed  by  the  FC,  battery  and  UC,  respectively.  To  sufficiently  prove  validity  of  the  experimental  setup,  the 
measured  UC  power  is  compared  to  the  calculated  UC  power  in  Fig.  26  with  the  same  scale.  High  similarity  between  the 
experimental  result  and  the  theoretical  value  exists  not  only  during  the  positive  power  transferring  process  (outputting  power)  but 
also  during  the  negative  power  transferring  process  (absorbing  power).  Consequently,  the  experimental  results  validate  that  the 
proposed  wavelet-transform-based  power  sharing  algorithm  is  applicable  to  the  real  hardware  environment. 
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Fig.  23.  Downscaled  total  power  demand  at  a  certain  time  interval  of  OCC.  Fig.  24.  Calculated  power  distribution  for  three  power  sources. 
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Fig.  25.  Experimental  results  of  power  distribution.  Fig.  26.  Comparison  of  UC  power  between  calculation  and  experimental  results. 


VII.  Conclusion 

With  the  proposed  wavelet-transform  based  power  management  algorithm,  the  decomposed  components  of  the  load  profile 
during  a  drive  cycle  can  be  distributed  to  the  FC,  battery  and  UC,  compatible  with  their  respective  characteristics  in  order  to 
obtain  satisfactory  performance.  The  lifetime  of  the  FC  module  can  be  significantly  increased  and  its  cost  can  be  dramatically 
reduced  by  using  the  architecture  and  algorithm  proposed  in  the  paper.  This  is  because  the  UC  deals  with  the  fast-variation 
transient  and  peak  power  demand,  while  the  battery  helps  the  FC  satisfy  the  low-frequency  portion,  with  the  latter  absorbing  the 
majority  of  the  power  demand  in  the  normal  running  process  except  acceleration  or  deceleration.  The  proposed  power 
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management  system  was  implemented  using  a  modeling  and  simulation  software  platform.  Experimental  data  has  been  validated 
against  the  simulated  data,  thus  verifying  the  accuracy  of  the  power  source  models.  The  simulation  results  also  verify  the  validity 
of  the  power-sharing  algorithm.  Real  time  controller  based  experiments  demonstrate  the  applicability  of  the  proposed  strategy  to 
the  real  system. 
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